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Abstract We present general relativistic radiation transfer formulations which 
include opacity effects due to absorption, emission and scattering explicitly. We 
consider a moment expansions for the transfer in the presence of scattering. The 
formulation is applied to calculation emissions from accretion and outflows in 
black-hole systems. Cases with thin accretion disks and accretion tori are consid- 
ered. Effects, such as emission anisotropy, non-stationary flows and geometrical 
self-occultation are investigated. Polarisation transfer in curved space-time is dis- 
cussed qualitatively. 
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1 INTRODUCTION 

Observations of the AGN MCG-6-30-15 by A SC A showed Fe Ka emission line with a very 
broad asymmetric profile (Tanaka et al. 1995). The presence of the line was confirmed by later 
observations using e.g. Chandra and XMM-Newton. The broad asymmetric Fe Ka lines have 
also been seen in the spectra of several other AGN (e.g. NGC3516, Nandra et al. 1999). The 
line is believed to be fluorescent emission from accretion flows very close to the black- hole event 
horizon, where gravitational and kinematic effects are important. There are several mechanisms 
contributing to the line broadening and profile modification. The rotational motion in the 
accretion disk causes the two peaks of the line; relativistic boosting makes the line asymmetric 
(with an enhanced blue (high-energy) line peak and a suppressed red (low-energy) peak); time 
dilation, due to gravity and transverse relativistic motion, broadens the line and shifts it to a 
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lower energy. (See e.g. review by Fabian et al. 2000.) Gravitational lensing also alters the line 
profile by affect the projection of the emission. 

Modeling relativistic lines from black-hole system have been carried out since the 70’s. In 
spite of the variations in the numerical/analytic techniques, line-profile calculations generally 
follow a prescription which can be summarised as follow. (1) Define the space-time metric. (2) 
Construct a model flat accretion disk to give a macroscopic velocity profile for the emitters. 
The disk can be stationary or time-dependent, Keplerian or non-Keplerian. Usually a para- 
metric emissivity profile is assumed. Local variations in the disk emission are modeled through 
parametric prescriptions of turbulence, or density /emissivity inhomogeneities, time-dependent 
flares etc. (3) Determine the geodesics of the “observable” photons and calculate the rel- 

ative energy shifts of the photons between the emitter and the observer. (4) Bin the photons 
according to their energies and sum the bins to produce the line spectrum. For time-dependent, 
line profile calculations, the photons are also time binned and summed, according to their ar- 
rival time. These calculations have provided insights for our understanding of relativistic effects 
on line profile in accreting black holes. However, as our observational capability improves, we 
need more sophisticate radiative transfer calculations for emission from black-hole systems in 
order to fully utilise the high quality spectral and timing data that will be collected by the next 
generation observatories. 


2 RADIATIVE TRANSFER FORMULATION 

Unless otherwise stated, we use the convention of natural unit (c = G = h = 1, where c is 
the speed of light, G is the gravitational constant, and h is the Planck constant). Without 
scattering, the general relativistic transfer equation reads 

= + (i) 

(Baschek et al. 1997; Fuerst & Wu 2004; Wu et al. 2006, see also Lindquist 1996), where A is the 
affine parameter, is the affine connection (determined by the metric g Ql 3 , which is specified 
by dr 2 = g a ^dx a dx (3 ) 1 k a is the 4-momentum of photons, and u a is the 4- velocity of the medium. 
The Lorentz invariant intensity 1 = 7„/i/ 3 , where I u is the specific intensity, v is the frequency, 
and \o and r) lt are the Lorentz invariant absorption and emission coefficients respectively. (The 
subscript u o” denotes that the variable is evaluated at the local rest frame.) This formulation 
allows the radiation to be self-absorbed internally in the emitting medium, and absorption 
and re-emission in external line-of-sight media. Once the photon geodesics are determined and 
the emission and absorption coefficients are given, the radiative- transfer equation can be solved 
along the line-of-sight (Fuerst & Wu 2004) using a ray-tracing, and the spectrum of the emission 
is determined. 

In the presence of scattering, photons are injected into the line-of-sight from a crossing ray. 
Equation (1) is therefore inadequate. In principle, adiative transfer in a scattering medium must 
be evaluated globally, instead of tracing individual rays independently each time. Obtaining 
a solution to the equation is, however, not straightforward and often practically impossible. 
Various techniques have been developed to bypass a direct evaluation of the proper integro- 
differential radiative transfer equation, for example, to obtain a solution to a set of differential 
equations, that approximate the exact radiative transfer equation. 

For general relativistic radiative transfer for emission from accreting black holes in AGN, 
the moment methods described in Fuerst (2005) is generally applicable. The method may be 
summarised as follow. First, expand the intensity Z(x Q , k a ) into a series using a set of orthogonal 
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symmetric tensors J Q , which gives 


l(x a ,k a ) = Oo J + £ai^ lQ 2 .... Q| n Q 'n a ’....n°' , (2) 

i«l 

where n a is the directional unit vector of the photon. Next, generate a set of projected tensors 
Jaia 2 .. .q» .... satisfying the condition: 


u Q J c 


orti«2 — Qt — 


oc 3a 


(3) 


The 4-velocity u Q is specified in a preferred reference frame, such as the local rest-frame of the 
medium. We then consider an approximation by truncating the moment expansion. For the j-th 
order, we have 


Xj(* Q ,* Q ) = , 


(4) 


where m a = n a + u Q . With these, we derive moment equations, whose of solutions 3aia 2 ...cn 
form the basis for the construction of various order approximations to the ‘true’ intensity 

For instance, if opacity of the medium is contributed by three processes: a free- free process 
due to a sub-population of thermal electrons (of weight X), synchrotron radiation due to another 
sub-population of electrons which is non-thermal (of weight (1 — X)), and electron scattering, 
the first-order moment equation will read 


m a 


= - \X<TH + (1 - X)cr syn + <r T | pj a m a + Xa tt pB + (1 - X)a syn pS liyn 


+ <J-TPJ a u a , 


(5) 


where cr tf.syn.T are the respective absorption cross section of the free-free, synchrotron and 
electron scattering processes, B = is the Lorentz invariant Planck function, S syn is the 

corresponding Lorentz invariant source function for synchrotron radiation, and the variable 

1 DF 

Z = ~p2 ~HX = “° n «;9^+ ,1 X/ (6) 


(with E = v for the photons). If the synchrotron emission is from a thermal population of 
electrons, the source function S is the same as the Planck function 3. The equation is analogous 
to the conventional radiative transfer equation. The first group of terms on the right side 
correspond to removable of photons by various processes; the next term corresponds to free-free 
emission; the term after it corresponds to synchrotron emission; the last term corresponds to 
injection of photons into the line-of-sight via scattering. The second-order and higher-order 
equations can be derived accordingly, following similar procedures (Fuerst 2005). 


3 THIN RELATIVISTIC ACCRETION DISKS 

The radiative transfer formulations presented in §2 is applicable to a wide range of astrophysical 
settings. For emission lines from accretion and outflows in black-hole systems, both formulations 
easily generate relativistic lines as those in Cunningham (1975), Fabian et. al. (1989), Stella 

(1990), Laor (1991), by setting the absorption, emission and scattering coefficient to zero. 

The inclusion of absorption, emission and scattering coefficients, the formulations take account 
of the full transfer effects. 

We show in the following subsection two examples of applications for geometrically thin 
relativistic accretion disks around black holes. The first one concerns angle-dependent effects 
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Fig. 1 The top panel of the left column shows t he comparison between a reflection spectrum 
of a geometrically thin relativistic accretion disk around a Schwarzschild black hole (dotted 
line) and a rest-frame reflection spectrum (solid line) generated by PEXRAV in XSPEC 
based on the model calculation of Magdziarz & Zdziarski (1995). The viewing inclination 
angle of the disk is 45°. The strength of the reflection follows a radial power-law profile with 
an index of —2. The bottom panel of the column shows the corresponding quotient spectra, 
w r hich is obtained by dividing the spectra by the rest-frame spectrum. The top panel of 
the right column show reflection spectra of thin accretion disks around black holes with spin 
parameters a = 0. 0.5 and 0.998 (solid, dotted and dot -dashed lines respectively). The bottom 
panel of the column shows the corresponding quotient spectra, in reference to the case with 
a = 0. In all the reflection spectra the intensity 1(E) is in an arbitrary unit. The quotient 
spectra are normalised, with the normalised intensity R = 1 at the energy E = 1 keV. 


on the reflection continuum. The second concerns a complex situation that the emitting and 
absorbing media vary with time and that emission and absorption are not confined to the disk 
equatorial plane, so that angle-dependent effects are treated explicitly throughout the entire 
region. 


3.1 angle-dependent emission - reflection spectra 

An AGN X-ray spectrum generally consists of a thermal black-body component, a hard power- 
law component, a reflection component, and lines, e.g. the fluorescent Fe Ka line. The thermal 
black-body component, the powder-law component and the lines Eire probably isotropic. The 
reflection component is, however, strongly viewing angle dependent. In the analyses of X-ray 
spectra data, model fits to the reflection component often assume a fix viewing angle, chosen 
to be the same as the viewing inclination of the normal of a geometrically thin accretion disk. 
This assumption can be problematic, as the space-time in the inner accretion disks is curved 
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and the emission is gravitational lensing. The pitch angles of the photons from a disk surface 
element that reach a distant observer are not the same as the inclination angle specifying the 
orientation of the accretion disk. Light aberration is also non- negligible as the flows in the inner 
accretion disk are highly relativistic. The situation is further complicated by the rotation of 
the black hole, which leads ^6 reference frame dragging. Moreover, limb effects occurs when the 
emission emerging from the disk surface. (See Lee et al. (2007) and Wu et al. (2006).) 

Gravitational and relativistic kinematic effects coupling with angle-dependence of reflection 
will lead to a number of observable consequences. Figure 1 illustrated some of the consequences 
by (i) a comparison between a relativistic disk reflection spectrum and a rest-frame spectrum, 
and (ii) a comparison between reflection spectra from accretion disks around black holes with 
different spin parameters. The first demonstrates how local angle-dependence together with 
relativistic energy shift modifying the reflection continuum. The second shows how the effects 
combine with the rotation of the black hole, giving rise to broad features mimicking relativistic 
emission lines. 

An obvious difference between a relativistic reflection spectrum and a rest-frame spectrum is 
the smearing of the absorption edges. The quotient spectra in Fig. 1 imply that fitting data of a 
reflection spectrum of a relativistic accretion disk using a model rest-frame reflection spectra will 
result in flux excess blue-ward of the edge and flux deficiency red-ward of the edge. These may 
lead to an artifact blue-shift emission line, which has a wavelength consistent with the element 
species associated with the edge. The artifact feature may also be mistaken as an emission line 
from relativistic outflows. Another obvious difference is at the 10 — 100 keV spectral region. 
The reflection peak is lowered and the high-energy tail is flattened in the relativistic spectrum. 
These will affect the slope of the fit power-law for the observed spectrum. The effect is severe 
in situations where the spectral coverage at the high-energies, say above 20 keV, is unavailable. 

The effects of the black hole’s rotation on the reflection spectra of relativistic disks are 
manifested predominantly in the red-ward region of the absorption edge (panels in the right 
column. Fig. 1. The disk reflection spectrum of a fast rotating black hole shows flux deficiency 
red-ward of the edge when comparing with the disk reflection spectrum of a slower rotating 
black hole. Using spectral models corresponding to fast rotating black hole in fits to the spectra 
data of a system with a slowly rotating black hole will therefore yield flux excess red- ward of the 
edge. The flux excess may create a line-like feature with a profile resembling those of relativistic 
emission lines for Schwarzschild black hole or even fast rotating Kerr black holes, depending 
the relative rotational rate of the two black holes in the comparison. The high-energy tail of 
the reflect spectrum is also affected by the black-hole spin. The general trend is that he more 
rapidly is the rotation of the black hole the steeper the slope of the high-energy spectral tail. 


3.2 time-varying accretion and outflow 

The current time-dependent calculations of emission from relativistic flows are mostly focused 
on effects of dynamics. They generally employ a generic ray-tracing algorithm which does not 
include opacity effects in the emitting and transmitting media. Although these studies have 
provided great insights on how space-time curvature and relativistic kinematics on X-ray spec- 
tral formation in accreting black holes, for proper calculations of emission from magnetised 
disk-outflow system, one needs to take full account of radiative transfer to treat the emission 
and absorption processes explicitly. 

The change of an emission spectrum with the viewing inclination of a system is not restricted 
to photon pitch-angle isotropy in reference to the normals of the emission surface elements. It 
also occurs when the emission or absorption of the line-of-sight material is anisotropic. For 
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Fig. 2 The panels on the top row show images of emissions from the flows, viewed at an 
inclination angle of 45°, at time </r* = 0 (left) and 300 (right) (where r, = r 8 /c, with r 8 is 
the Schwarzschild radius of the black hole), corresponding to an initial state and a later state 
of the simulation. Bright synchrotron filaments (in greenish blue colour) are seen near the 
disk surface. The bright blue ring at the center is first-order lensed disk image. The fuzzy 
yellowish ring near the center in the image at f/r 8 = 300 is due to direct thermal free-free 
emission from dense material piled up in the inner accretion disk. The panels on the bottom 
row show images of emissions from flows viewed at an inclination angle of 85° for t/r a = 0 
(left) and 300 (right). The faint vertical features below the disk are the components of the 
counterpart jet. 


instance, synchrotron radiation is strongly beamed intrinsically, usually in the direction along 
the magnetic field line. 

Magnetised accretion disks and outflows are known to be non-steady. The magnetic fields 
threading the flows are twisted continually and the local field strengths fluctuating. The rapid 
variations in the field configuration and in the flow coupling with the gravitational lensing give 
rise to a number of unexpected phenomena. Figure 2 shows emission from from a flat accretion 
disk and its outflow. The dynamical variables are obtained by 2.5-D general relativistic MHD 
simulations using the 3-D RAISHIN code (Mizuno et al. 2006a,b). In the radiative transfer 
calculations the opacity is assumed to be due to thermal free-free and synchrotron processes. 
Bright filaments are seen in the disk image. They are caused by chance alignment of magnetic 
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Fig. 3 Energy shift images of opaque relativistic accretion tori around a rotating black 
hole with a spin parameter a = 0.998, viewing at inclination angles i = 1°, 30°, 60° and 
89° (panels ordered from left to right and top to bottom). Red for red shifts in energy and 
blue for blue shifts. The tori have a rotational velocity profile index n = 0.1 (see Fuerst 
&: Wu 2007) and the innermost boundary of the emission surface reaches 1.3r g , where the 
gravitation radius r g = 1 for normalised black-hole mass. 


field lines, along which synchrotron emission is beamed and hence boosted. As the simulations 
evolve, the filaments wiggle quasi-periodically (i.e. QPOs) on time scales of the order that of 
the dynamics of local flows. The filament features are robust — seen at all viewing inclinations. 
They are geometrical in nature, and are not sheared by the differential rotation of the accretion 
disk, This is in contrast to material clumps circulating around in the accretion disk, which 
disperse after a few orbits. Note that there is a pair of jets, with the counter-jet below the 
equatorial plane is visible. The simulations show two concentric components in each jet, each 
of which have different dynamical and emission properties, and a material piled up in the inner 
disk, giving rise to strong free-free emission. 
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Fig. 4 The left panel shows the emission line profiles of an opaque torus with parameters as 
those in Fig. 3 (i.e. n = 0.100). viewed at inclination angles of i = 1°, 15° , 30°, 45°, 60°. 75° 
and 89°. The right panel shows the emission line profiles of an opaque torus with parameters 
the same as those of the torus in the left panel except that n = 0.001 (i.e. practically a thin 
disk). The intensity 1(E) is normalised such that 1(E) = 1 at E/E Q = 1 for » = 45°. 


4 RELATIVISTIC ACCRETION TORI 

4.1 Opaque tori 

An accretion torus is a 3-D object, unlike a geometrically thin accretion disk, which is 2-D. 
Several methods commonly used in relativistic line calculations for thin accretion disks around 
black holes (e.g. Cunningham 1975) are not applicable, as for tori the emission originates from 
emission surfaces not in the equatorial plane. The situation is further complicated if the emission 
is angle dependent. While the entire upper surface of a geometrical thin accretion disk is always 
visible, self-occ illation effects in tori are important at high viewing inclination angles (Fig. 3). 
This plays a very significant role in shaping the emission line profiles. Figure 4 shows the profiles 
of emission lines from two accretion tori with different thickness (specified by the velocity profile 
parameter, n, Fuerst & Wu 2007) viewing at various inclination angles. The torus with n = 0.1 
has substantial thickness (see Fig. 3); the torus with n = 0.001 is practically the same as a 
geometrically thin disk. Both tori show similar line profiles at low inclination angles . The lines 
are red-shifted because of gravitational red-shift and transverse Doppler effect. The difference 
between lines from the two tori become more visible for viewing inclinations i larger than ~ 45°. 
At high inclination angles, the lines from tori with the larger n are narrower, i.e. suffer less red- 
shuft and less blue shift, as the innermost region, where relativistic effects are more severe, is not 
self-obscured. Here, it demonstrates that extracting information of the parameter of the system 
is non-trivial, as geometrical factors (such as aspect ratio of the accretion torus/disk) could 
affect maximum red-shift and maximum blue-shift in the line profile, thus causing confusing in 
the estimation of the location large stable particle orbit and hence the spin of the black hole. 

4.2 Semi-transparent tori 

Another different between an accretion torus and an accretion disk is that the emission is 
not restricted to a surface when the torus is not optically thin. In this situation, the emis- 
sion seen by the distant observer is contributed by the entire accretion torus, as shown in 
Figure. 5). The relativistic beaming is more important than for that in opaque tori, where the 
most beamed/ boosted emission is hidden because of self-occulation. If the torus is semi-opaque, 
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Fig. 5 Surface brightness of semi-transparent accretion tori around black holes with a = 
0.998, viewed at inclinations i = 15°, 45°, 60° and 85° (panels from left to right amd from top 
to bottom). The opacity is provided by the Fe K a and K (3 lines. The brightness is normalised 
such that the maximum pixel brightnesses of all tori in the images are the same. 


various complicated opacity effects may arise. For example, two lines with small difference in 
wavelengths, the line emitted from one part of the torus may be resonant with another line at 
the other part of the torus, because of the energy shifts during to velocity shears and difference 
in the gravitational potential of the emitters and the line-of-sight absorber/scatterer. While the 
emission lines from such tori are still broad and asymmetric as those of the opaque disks and 
torus, the line profiles are not at identical (Fuerst 2005). 


4.3 Scattering dominated tori 

In the scattering torus, the opacity is provided by electron scattering as well as the 
line/continuum emission and absorption processes. As photons can be scattered into a ray 
from a cross ray, the radiative transfer needs to be evaluated globally. With the moment for- 
mulation described in §2, the emission from scattering accretion tori can be calculated, and 
Figure 6 shows an example. Although a scattering torus may resemble a semi-transparent ab- 
sorptive torus, the emission properties are very different. Firstly, line resonance is difficult to 
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Fig. 6 Surface brightness of scattering accretion tori around black holes with a = 0.998, 
viewed at inclinations i = 45° (left) and 85° (right). The opacity is due to eleectron scattering. 
The brightness is normalised such that the maximum pixel brightnesses of the two tori are 
the same. 


occur for a strongly scattering torus. Secondly, the beaming effects less obvious in a torus with 
electron scattering, as the beamed emission can be scattered out from the line-of-sight. This is 
in big contrast to to an absorptive torus, where a beamed line photon is highly energy-shifted 
that it will not suffer being resonant scattered/absorbed. The total emission from an opaque 
torus is the sum of weighted and shifted Planck function from the visible torus surface element. 
This is not true for a thick torus, the emission from a scattering thick torus is not necessarily 
thermal and the total emission from it is not simply summing the weighted and shifted local 
Planck function of the visible elements of the last scattering torus surface. In addition to all 
these, an opaque absorptive torus would not show strong polarisation, but a scattering thick 
torus is expected to be polarised. 


5 BRIEF COMMENTS ON POLARISATION RADIATIVE TRANSFER 

In Newtonian space-time, the polarised radiative transfer equation may be expressed as 

i> /,M) = -aVjM) +*(«', fi) + / d "'/ . (?) 

where i and j are indices for the polarisation modes, l\ is the intensity of a polarisation mode, 
ji is the emission coefficient, K §j is the opacity whose diagonal elements specify absorption of 
the polarisation and off-diagonal elements specific polarisation conversion and rotation, Sjj is 
the scattering matrix, v is the frequency of the radiation, and Cl is the directional unit vector 
of the radiation propagation. The propagation differential operator is given by 

V = iJUn.v. 
cat 


Equation (7) implies that the creation of a polariation mode and the conversion between differ- 
ent polarisation modes are caused by the meida, either the emitters or the line-of-sight material. 
For instance, in the presence of a magnetic field, gyrating electrons will give rise to cyclotron and 
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“torsion" along the propagation 



Fig. 7 Schematic illustrations showing the rotation of polarisation vectors in the presence 
of magnetic field (left) and in rotational space-time (right) as the radiation propagates. 


synchrotron radiations, which are polarised. Also, the polarisation will undergo Faraday rota- 
tion/conversion when the radiations propagate in the magnetised plasma (see Fig. 7). Scattering 
processes may also induce polarisation. The propagation operator, however, does not affect the 
polarisation properties of the radiation. In curved space-time, this may not be true. 3-vectors 
defined on a local Lorentzian space-time are not preserved undergoing parallel transport along 
a geodesic (see e.g. Misner, Thorne &: Wheeler 1973). As the radiation propagates in a curved 
space-time, the orientations and strengths of the electric and magnetic fields change as seen by 
a distant observer, resulting in rotation, stretching and twisting of the polarisation vectors. If 
the space-time is rotating, such as that around a Kerr black hole, the polarisation vector can 
also be dragged into rotation (Fig. 7). For radiative transfer in black-hole environments, con- 
version and rotation of the polarisation vector are caused by space-time curvature and rotation 
as well as the radiative properties of the emitters and the media through which the radiation 
propagates. To disentangle all these effects are non-trivial. However, before this difficult task 
is carried out, one needs a fully covariant polarisation radiative transfer formulation applicable 
and suitable for astrophysical settings, which is currently lacking. 

6 SUMMARY 

We present a general relativistic formulations for radiative transfer which treats emission and 
absorption processes explicitly. We also present a more general formulation which includes scat- 
tering in addition to emission and absorption. A moment solution is derived for this equation. 
The formulations are applicable for the calculations of emission from accretion disks and accre- 
tion tori around black holes. Cases of reflective flat disks, opaque tori, semi-transparent tori and 
scattering tori are presented to illustrate various effects, which have been failed to recognise by 
conventional ray-tracing relativistic line calculations. Polarisation transfer in curved space-time 
is briefly discussed. 

We summary our findings as follow. (1) Angle-dependent effects are important and im- 
proper modeling of the reflection continuum will give rise to artifact features with profiles very 
similar to those of relativistic lines from accretion disks around fast rotating black holes. (2) 
Time dependent flows in the presence of magnetic fields together with relativistic beaming can 
give rise to QPOs in the emission. (3) Self-occulation is important for opaque accretion disks 
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with substantial thickness (i.e. accretion tori) at high viewing inclination angles. Lines from the 
thicker tori are narrower and show less relativistic effects, as emission from the inner surface, 
where the gravitational field is the strongest and the flow is the faster, is obscured. (3) Emission 
from semi-transparent accretion tori show more relativistic beaming effects. Processes forbid- 
den in opaque disks/tori, sucL as resonance of lines with different wavelengths, can occur in 
semi-transparent disks, due to gravitational potential gradient and velocity shear. (4) Opaque 
scattering tori and opaque absorptive tori have different radiative properties, and an opaque 
scattering torus may not have a differential, weighted and relativistic energy-shifted thermal 
spectrum. 
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